arXiv:astro-ph/9706217vl 20 Jun 1997 


The Relation Between Fourier Coefficients and RR Lyrae Magnitudes 

Ann Berit Saust 

Princeton University Observatory 
Princeton, NJ 08544-1001, USA. 

ABSTRACT 

The relation between the Fourier coefficients determined by the light curves of RR Lyrae 
variables in the Sculptor dwarf galaxy, the M5 cluster, and the unique globular cluster uj 
Centauri is investigated. A couple of recent papers claim that it is possible to determine both 
metallicity and absolute magnitude of RR Lyrae stars from their period and two of the Fourier 
coefficients. However, this investigation show that fitting high order Fourier series to noisy 
and/or sparse data results in a large scatter in the determination of the absolute magnitudes. 
Unless we can find a reasonable way of smoothing noisy data, it will be difficult to use RR 
Lyrae’s as standard candles. 


Subject headings: stars: variables: other (RR Lyrae) — stars: fundamental parameters — 
globular clusters: individual (M5, uj Centauri) — galaxies: individual (Sculptor) 


1. Introduction 

It has long been believed that the absolute magnitude of RR Lyrae stars could be determined by their 
light curves and, to some still unknown extent, their metallicity. Fourier decompositions, were introduced 
by Simon & Lee 1981, as a means of qualitative description of light curves of pulsation variables in general, 
and they showed that amplitude ratios and phase differences provide a useful description of the Hertzsprung 
progression for classical Cepheids. In addition, Simon & Teays 1982 showed that the Fourier decomposition 
parameters of 70 RR Lyrae field stars are more sensitive discriminators of the Bailey type (ab or c) than 
the traditionally employed period-amplitude diagram. 

Several globular clusters (e.g. uj Centauri) contain many RR Lyrae stars, and since the variables in a 
globular cluster probably constitute a much more uniform sample than the field RR Lyrae’s, a considerably 
smaller scatter is expected in data from one cluster. Petersen 1984 examined Fourier decompositions of RR 
Lyrae’s in uj Centauri, and found that the scatter in the Fourier parameters was larger than expected for 
a uniform sample. Petersen 1984 concluded that the scatter was most likely due to the large metallicity 
distribution ([FE/H] ranges from -2.3 to -0.5 (Butler et al. 1978), differences in mass of the horizontal 
branch stars, or differences in effective temperature. As Smith 1995 pointed out, uj Centauri is so unusual 
in its range of chemical composition that it is questionable whether uj Centauri is the key to the absolute 
magnitude-metallicity relation, too unusual to be a representative of the RR Lyrae population, or just 
one more clue to an absolute magnitude-metallicity relation relationship which is more complicated than 
expected. See Smith 1995 for a review on RR Lyrae stars in general. 

According to two recent papers, Kovacs & Jurcsik 1996 (KJ) and Jurcsik & Kovacs 1996 (JK), the 
metallicity and absolute magnitude of RR Lyrae ab stars can be described by a simple linear combination 
of the period and two Fourier coefficients. Kovacs & Jurcsik 1997 added multi band observations to their 
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earlier results to obtain an even smaller error, leading to an estimate of the relative distance moduli with 
an accuracy of < 0.03 mag. However, their sample did not include the unique globular cluster oj Centauri, 
which is included here. In addition to to Centauri, this investigation includes the Sculptor galaxy and the 
M5 cluster which were also included in JK and KJ, such that results and errors can be compared directly. 


2. The Observations 

The observations of Sculptor galaxy and ui Centauri were made by the OGLE team during 1992 and 
1993 using the 1-m Swope telescope at Las Campanas Observatory (see Kaluzny et al. 1995 and Kaluzny 
et al. 1997). The OGLE data is available via anonymous FTP from astro.princeton.edu, in the directories 
bp/Sculptor and bp/Omega_Cen. All The M5 data were made by N. Reid during 1991 and 1992 using the 
the 60 inch telescope at Palomar Observatory (see Reid 1996) and is also available via anonymous FTP 
from astro.princeton.edu, in the directory bp/M5. All observations used in this paper were done in Johnson 
V. 


3. Fourier Analysis 


It is customary to decompose the light curve into its Nth order Fourier decomposition: 

N 

My = ao + afesm(27rfc/) + bkCos{2'Kkf) (1) 

fe=i 

where f is the phase given by / = {JDobs — JDo)/P, JDq being the mean epoch and P the period. 
Defining + b^, tan{(j)k) = — ak/bk this can be written: 

N 

My = ao + '^Hksin{27rkf + (j>k) (2) 

fe=i 

(see Petersen 1986 for more details). 

Since the high noise level in especially the OGLE data makes it impossible to make a direct high order 
Fourier fit, all the data was smoothed using polynomials. Since the same mechanism is believed to be 
causing the variability in RR Lyrae’s, it would be reasonable that the same smoothing technique could be 
applied to all the data. This also has the advantage that the results are reproducible and can be compared 
to future observations. 


The data were first wrapped according to the period, divided into 9 equal parts which were then 
smoothed with a 7th order polynomial. The high order of the polynomial ensures that the curve is smooth 
in intervals with many data points, but follows the data points in less crowded intervals. To ensure 
continuity, the actual fitting was done over a three times larger interval. This combination of windowing 
and order of the polynomial gave the smallest for all stars, where 
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where Nobs is the number of observations, iJ,obs,i is the i’th observed magnitude, is the i’th calculated 
magnitude, and ai is the noise. A high (8th) order Fourier series was then fitted to the smoothed light 
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curve using singular value decomposition as described in in Press et al. (1992). Just as in KJ, the Fourier 
parameters refer to a sine decomposition and the phase is chosen as the closest value to 5.1. Increasing the 
order of the Fourier series did not show any significant changes of the (lower) Fourier coefficients. 

Fig. 0 shows the observations and the Fourier fit for one star from the Sculptor galaxy. As can be seen, 
a better fit could be obtained with human aid, e.g. using different order of the polynomials in each interval, 
changing the size of the intervals, etc. (which is what KJ and JK did in their papers), but as mentioned 
above, a general smoothing technique has many advantages. No attempts were made to delete or disregard 
obvious bad data points on the individual light curves since this would be impossible in intervals with 
sparse data, and one might worry about introducing non random bias into the data. As a result, the scatter 
is expected to be fairly large; however, given sufficient amounts of data the final result should be similar to 
KJ’s. 

JK’s main conclusion is that the Fourier parameter defined as (^^31 = — 3 • (in general, 

4>ji = ~ J ■ '(’ 1)1 is the most important parameter when it comes determining the metallicity and KJ 

found that the three parameters, the period, P, iJi, and 4>si, completely determines both magnitude and 
metallicity. Their formula for the absolute magnitude is: 

Mv = 1.221 - 1.396P - 0.477iJi + 0.103(()3i (4) 

and the metallicity: 

[Fe/H] = - 5.038 - 5.394P + 1.345031 (5) 


Following KJ, the observed magnitudes were fitted as a linear function of P, iJi, and 03i, again using 
singular value decomposition. To avoid zero point calculations and problems, the three clusters are treated 
separately, and since RR Lyrae ab stars are both the most abundant in the sample and the only stars 
treated by JK and KJ, these stars are the only ones included in the following. In addition, stars showing 
peculiar behavior are omitted. This includes stars showing Blazhko behavior, stars with undetermined 
periods, a few stars with periods longer than 1 day, and stars several magnitudes fainter than the majority 
of the cluster. This means that 124 stars in the Sculptor galaxy, 35 M5 stars, and 87 out of 97 oj Centauri 
RR Lyrae ab stars were used. 

To estimate how good the fit is, is calculated for each ht, where here is defined as: 
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where Ngtar is the number of stars, mobs,i is the i’th observed mean magnitude, is the i’th calculated 

magnitude, and amean.i is the 1 sigma error on the observed mean magnitude. 

Fitting for the three parameters gives = 800 for Sculptor, = 77 for M5, and = 750 for oj 
Centauri. Fig. || shows the fits which are given by 


ms = 20.45 - 0.62P + 2.96 • lO-^Pi + 7.26-10-3031 (7) 

WM = 15.94 - 1.29P - 7.98-10-3Pi + 1.47-10-3031 (8) 

mo = 14.49 - 1.29 - IQ-^P + 7.06 - lO-^Pi - 9.33 - 10-3031 (9) 


As can be seen both from the figures and the x^ ’s, the hts are poor, much poorer than what was found 
in KJ. Adding additional Fourier coefficients to the fit gives only a marginal improvement. Since some of 
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the Fourier fits show large these were omitted in an attempt to reduce the scatter. However, neither 
the scatter nor the fit changed significantly, implying that the main contribution to the scatter is coming 
from bad data points in intervals with sparse data (assuming that it is possible to determine the magnitude 
as a function of P, Hi, and (j) 3 i in the first place). Notice also that there is almost no dependency on ^ 31 , 
in contrast to KJ’s result, and that the dependency on Hi is strongly influenced by the scatter. It should 
also be noted that while the fits for the Sculptor galaxy and the M5 cluster are somewhat similar, the fit for 
uj Centauri deviates significantly. The noise level in the two OGLE data sets are much larger than the noise 
level in most of KJ’s data, but even in the case of M5 which has a much lower noise, the ^31 parameter 
seems to be negligible. 

Different smoothing techniques (i.e. different windowing or orders of the fitted polynomials) did not 
change the fitted parameters significantly (deviations within 20%). Since the noise level in the M5 data is 
low, the M5 data was recalculated without any smoothing. Fig. ^ shows the result, is 93 here, and the 
fit is given by: 

niM = 15.79 - 1.09P - l.ie-lO-^iJi - 6.14 • 10“3(()3i (10) 

Again, the dependence on the ^31 parameter is small. 

Since the metallicity distribution in uj Centauri is larger than in any other cluster, it would be 
interesting to do the same exercise for the metallicity. Unfortunately, there is not enough data available in 
the literature, and the metallicity of the to Centauri stars were instead calculated using the result from JK 
(equation || above) and compared to the observed [Fe/H] values from Butler et al. 1978. Notice that there 
are only 18 stars which are common for the two data sets. The result is shown in fig. ^ and, as can be seen, 
resembles a scatter plot. 


4. Discussion 

Although the observations in this paper are noisy, the deviations from earlier results are troublesome 
in both the case of Sculptor and M5, where it was impossible to reproduce the low scatter found by KJ, 
and in the case of lu Centauri where none of the earlier fits seems to fit the observations. 

Since KJ also included Sculptor and M5 in their analysis, their fits should show scatter similar to what 
was found above. It seems that the deviations for Sculptor and M5 are mostly caused by the different 
smoothing techniques. However, any general smoothing using polynomials all give results similar to the 
above, and looking at the individual light curves, any smoothing (general or done with human aid) is a non 
trivial task due to the noisy and/or sparse data. In addition, the fits for both Sculptor and uj Centauri (fig. 
H) shows evidence of extinction. 

Due to the unique metallicity distribution in lo Centauri this cluster can be regarded as either the best 
cluster to determine relations between absolute magnitude and metallicity or a non representative example 
of the RR Lyrae population. This investigation seems to support the latter; however, the concerns over the 
smoothing makes it impossible to make a conclusion. 

Since noisy and sparse data is common place, a general smoothing technique would be desirable, and 
it seems that if we want to use RR Lyrae’s as standard candles, some form of general consensus on the 
smoothing has to be reached. It is by no means certain that the smoothing applied here is the best available, 
and it would not be surprising if the scatter is reduced by a different way of smoothing. However, based 
on the above, it seems unlikely that any general smoothing would reduce the scatter to the level found by 
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Fig. 3.— Same as fig. 2 above for M5, only no smoothing was applied. 



Fig. 4.— Observed and calculated [Fe/H] for oj Centauri. [Fe/H] was calculated according to the formula 
given by Jurcsik & Kovacs 1996. 
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KJ. Since the above results strongly suggest that the main part of the scatter is due to bad data points in 
intervals with sparse data, a better sampling of data would obviously improve the results considerably. 
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